9, 2009 13:44 WSPC/INSTRUCTION FILE ijqianncaling'r 



International Journal of Quantum Information 
© World Scientific Publishing Company 



NON-HERMITIAN ADIABATIC QUANTUM OPTIMIZATION 



GENNADY P. BERMAN 

Theoretical Division, MS B213, Los Alamos National Laboratory 
Los Alamos, New Mexico 87544, USA 
gpb@lanl.gov 

ALEXANDER I. NESTEROV 

Departamento de Fisica, CUCEI, Universidad de Guadalajara, Av. Revolucion 1500 
Guadalajara, Jalisco, 444&0, Mexico 
nesterov&cencar. udg.mx 

Received Day Month Year 
Revised Day Month Year 

We propose a novel non-Hermitian adiabatic quantum optimization algorithm. One of 
the new ideas is to use a non-Hermitian auxiliary "initial" Hamiltonian that provides an 
effective level repulsion for the main Hamiltonian. This effect enables us to develop an 
adiabatic theory which determines ground state much more efficiently than Hcrmitian 
methods. 
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1. Introduction 

Many physical and combinatorial problems associated with complex networks of 
interacting degrees of freedom can be mapped to equivalent problems of finding 
the ground (or minimum cost) state of a corresponding quantum Hamiltonian Ho 
| 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | Q ne Q f j- ne approaches to finding the ground state of Ho is adia- 
batic quantum computation which can be formulated as follows. Consider the time 
dependent Hamiltonian H(t) — (t/T)TCo + (1 — t/r)Hx, where Ho is the Hamilto- 
nian whose ground state is to be found, r is the given time-interval of quantum 
computation, Hi is an auxiliary "initial" Hamiltonian and [7io,7ii] ^ 0. As time 
varies from t = to t = r, the Hamiltonian interpolates between Hi and Ho- 

If the system is initially close to the ground state of Hi, and if r is sufficiently 
large (slow variation) , then the system will remain close to the instantaneous ground 
state (i.e., that of H T (t)) for all t £ [0, r]. In particular, at t — r the ground state 
of the total Hamiltonian, H T , will be close to the ground state of Ho, which is the 
state we seek. In practice, Hi is chosen such that its ground state is known, then 
the dynamics is allowed to evolve and the state of the system evolves into the final 
state which is the solution to the problem. 
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Closely related to the adiabatic quantum computation is quantum annealing, in 
which one has a Hamiltonian Ho to be optimized, and an auxiliary (non-commuting) 
term Hi so that the total Hamiltonian reads Htot{t) — Ho + F(t)Hi, where T(t) 
reduces from very high value to zero during the evolution. The coefficient T(t) is 
the control parameter and initially T kept very high so that Hi dominates over Ho- 
One starts with the ground state of Hi as the initial state, and if T(t) is slowly 
decreasing, the adiabatic theorem guarantees reaching of the ground state of the 
main Hamiltonian Ho at the end of computation, assuming that there are no energy 
level crossings between the ground and excited states. 

Quantum optimization algorithms usually require the presence of a gap between 
the ground state and first excited state, however, in typical cases the minimal gap, 
g m , is exponentially small 0. This increases drastically the total computational time 
and from a practical point of view, the advantage of the method is lost. 

In this paper, we propose a novel adiabatic non-Hermitian quantum optimiza- 
tion algorithm. The non-Hermitian Hamiltonians naturally appear when a quan- 
tum system has both discrete (intrinsic) and continuous spectra, and one per- 
forms a proiection of the total wave function on the discrete part of the spec- 
trum | 10 | 11 | 12 | 13 | 14 | 15 |16| j n casej corresponding intrinsic energy levels 
acquire the finite widths, which are associated with the transitions from the in- 
trinsic states to the continuum. Then the dynamics of the intrinsic states can be 
desc ribed by the S chrodinger equation with an effective non-Hermitian Hamilto- 
nian EUSHZHSHSl The adequate candidates for implementation of non-Hermitian 
adiabatic qua ntum optimiza tion algorithm can be, for example, superconducting 
phase quMtsBS MESmSl 

We show that coupling the system to a non-Hermitian auxiliary initial Hamil- 
tonian induces an effective level repulsion for the total Hamiltonian. This effect 
enables us to develop an adiabatic theory without the usual gap restriction, and to 
determine much more efficiently the ground state of Ho- 

2. Adiabatic Quantum Optimization 

The generic adiabatic quantum optimization problem may be formulated as follows 

3. Let Ho be a Hamiltonian whose ground state is to be found, and Hi be the 
auxiliary "initial" Hamiltonian. Then we consider the following time -dependent 
Hamiltonian: 

H T (t)=fo(t)Ho + fl(t)Hl, (1) 

where {Ho,Hi] =/= . The functions fo(t) and fi(t) are monotonic increasing and 
decreasing, respectively, and satisfy the following conditions: /i(r) = 0, /o(t) = 1, 
and /i(0) > / (0), if / (0) / 0. The Hamiltonian H T {t) -> H Q , as t -> r, and we 

a For instance, in the commonly us ed quantum optimization n-qubit models, the estimate of the 
minimal energy gap is g m oc 2~ n / 2 H| 4 | 38 | 39 |40] 
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assume that, for any choice of the function fo(t), TL T {t) is dominated by TL\ at the 
initial time t = . 

The evolution of the system is determined by the Schrodinger equation 

ij^(*)>=«r(*M*)>- (2) 

We impose the initial conditions as follows: H.i\ijj g ) — E g \ip g ), where E g denotes 
the energy of the ground state \ip g ) = |?/>(0)), which is assumed to be the ground 
state of the auxiliary Hamiltonian Ti\. The adiabatic theorem guarantees that the 
initial state \ip g ) evolves into the final state \ip g (T)), which is the ground state of 
the Hamiltonian Hq as long as the instantaneous ground state of TC T (t) does not 
become degenerate at any time. 

The valid ity of the adiabatic theorem requires that the following condition be 
satisfied EE 



E 



(E m (t) - E n {t)f 



« 1, (3) 



where \ip n (t)) is the ground instantaneous state related to the instantaneous 
eigenenegry E n (t). 

The condition ^ is violated near the degeneracies in which the eigenvalues co- 
alesce. In most common case of double degeneracy with two linearly independent 
eigenvectors, the energy surfaces form the sheets of a double cone. The apex of 
the cones is called a "diabolic point" , and since, for a generic Hermitian Hamil- 
tonian, the co-di mensio n of the diabolic point is three, it can be characterized by 
three parameters | 26 | 27 ] Moreover, in the vicinity of the degeneracy point , th e N- 
dimensional problem effectively becomes a two-dimensional problem This 
will be essential in the following argument. 

For quantu m optim ization the commonly used version of the adiabatic theorem 
takes the forml 3 ! 4 ! 31 ! 

max\(MtWAtM g (t))\ U) 
mm\E e (t) - E g (t)\z ' [) 

where \^p g {t)) and \ip e (t)) are the ground instantaneous state and the first excited 
state of the total system, and "dot" denotes the derivative with respect to the 
dimensionless time s — t/r, TL T = dTL T /ds = TdH T (t)/dt. 

As can be observed, if during the evolution the gap AE(t) = \E e (t) — E g (t)\ 
becomes small enough, the amount of time required to pass from the initial state to 
the final state becomes very large and from the practical point of view the adiabatic 
quantum optimization loses its advantage compared with thermal annealing. 

Since in the vicinity of the level crossing point only the two-dimensional Jordan 
block related to the level crossing makes the most considerable contribution to the 
quantum evolution, the iV-dimensional problem can be describ ed by an effective 
two-dimensional Hamiltonian which can be obtained as follows S^l Let t c be the 
crossover point at which the energy gap between the ground state and the first 
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excited state of the total Hamiltonian TC T (t) achieves its minimum. In the two- 
dimensional subspace corresponding to E g (t c ) and E e (t c ), we choose an orthonor- 
mal basis {|0), |1)} and complement it to the complete basis of the iV-dimensional 
Hilbert space by adding the eigenvectors \ipk{t c )) {k — 2, . . . , N — 1). 
Now, an arbitrary state \ip(t)) can be expanded as 

N-i 

\m) = a(t)\0) + 0(t)\l) + J2 Ck(t)\MQ). (5) 

k 5^0,1 

Inserting this expansion into the Schrodinger equation @, we obtain the coeffi- 
cients, a(t) and 0(t), as the solution of the two-dimensional Schrodinger equation 

*^|u(t)>=tt e /(t)|u(t)>, (6) 

where 

nef[t) -{x(t) +l Y(t) X(t)-Z(t)J (7) 
and \u(t)) — ( a }f\ J • The matrix elements in Eq. ([7]) are determined by 



A(t) = i((0|W T (t)|0} + (l|W T (t)|l)), (8) 
X(t) = i((0|W T (t)|l} + (l|W T (t)|0)), (9) 

^) = i((0|W r (t)|0}-(l|W T (t)|l)). (11) 

Solving the characteristic equation for H e f(t), we obtain (below, we do not 
indicate in some expressions the explicit dependencies on t) 

E± = A ± VX 2 + Y 2 + Z 2 . (12) 

Setting R = (X, Y, Z) , we find the energy gap between the ground state and the 
first excited state is AE = 2R. 

Inserting H T (t) = f (i)Ho + fi(t)Hi into Eqs. (8j) - fTTj). wc can write the 
effective Hamiltonian as 

H ef (t) = A(t)l + / (t)-Ro • <t + h(t)Ri ■ a, (13) 

where \{t) = f (t)X + f 1 (t)X 1 , Rv = {X 0} Y ,Z„) and R l = [X X ,Y X ,Z X ). The 
time- independent parameters Aq, Ai and components of the vectors Rq and R x arc 
determined from Eqs. (jSJ) - PT|) by substitution of Ho and Hi instead of H T . Next, 
setting J = fo(t)Ro and g = fx(t)Ri, we obtain R = \/ g 2 — 2gJcosa + J 2 , where 
cos a = — Rq ■ Ri/RqRi. 
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At the crossover point we obtain 

j c - 5r c cosa 

g c = -J c - t , (14) 

g c — J c cos a 

where we denote g c = g(t c ), J c = J{t c ), g c = g(t c ) and J c = J(t c ). This yields 



' gl - 2g c J c cos a + J% 

|A£| min = -* ■ , ■ 2\J C \ sina. (15) 

\g c - J c cos a| 

It follows that sina ~ g m /2\J c \, where g m — |A_E| m i n is the minimum gap between 
the first two energy levels of the total Hamiltonian H T {t). 



3. Non-Hermitian Adiabatic Quantum Optimization 

Non-Hermitian quantum optimization can be implemented by the following gener- 
alization of the Hermitian adiabatic quantum optimization: 

Hr(t) = /o(t)«o + /i(*)fti, (16) 

where fi(t) — fi(t)—if2(t). The functions fo(t) and fi(t) (being monotonic decreas- 
ing and increasing, respectively) satisfy the following conditions: /i(0) ^> /o(0), if 
/o(0) 7^ 0, and /i(r) = 0, /o(t) = 1. In addition, we assume that the function f2(t) 
is monotonic and /2(t) = 0. 

The evolution of the to tal system is determined by the Schrodinger equation 
and its adjoint equation ^31 

i^ t m))=n r {t)m)), (17) 
-i^m)\ = m)\H T {t). (is) 

We impose the initial conditions as follows: ?ti\ipg) = E g \ip g ), E g being the energy 
of the ground state of the initial non-Hermitian Hamiltonian Hi = fi(0)H%. 

We denote by \ipn{t)) and (ip n (t)\ the right/left instantaneous eigenvectors of 
the total Hamiltonian: 

Hr(t)\Mt)) = En(t)\ll> n (t)), (19) 

$ n (t)\H T (t) = $ n (t)\E n (t). (20) 

We assume that both systems of left and right eigenvectors form a bi-orthonormal 
basis, {i> m {t)\i> n {t)) =5 mn ^. 

For the non-Hermitian quantum optimizationproblem the criterion validity of 
adiabatic approximation can be written ag | 33 | 35 |56] 



max\(^ e (t)\n r (t)\ip g (t))\ 



T>> mm\E e (t) - E g {t)? ' ^ 
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This restriction is violated near the ground state degeneracy where the complex 
energy levels cross. The point of degeneracy is known as the exceptional point, 
and it is characterized by a coalescence of eigenvalues and their corresponding 
eigenvectors, as well. Therefore, studying the behavior of the system in the vicinity 
of the exceptional point requires special care I 29 | 30 | 37 [ 

At the crossover point, t c , we introduce the bi-orthonormal basis as follows: In 
the two-dimensional subspace spanned by the ground state and the first excited 
state of the total non-Hermitian Hamiltonian, H T , we choose an orthonormal ba- 
sis {|0), |1)} and complement it to the complete basis by adding the eigenvectors 
\MQ) {k = 2,..., N-l). Then the set of states {(0|, (1|, $ k (t c )\; |0), |1), |Vfc(*o)», 
where k runs from 2 to N — 1, forms the bi-orthonormal basis of the iV-dimensional 
Hilbert space. 

Basically repeating the same procedure as for the Hermitian case, we obtain the 
effective two-dimensional non-Hermitian Hamiltonian as 

H ef (t) = A(t)l + f (t)Ro • o- + h{t)Ri ■ <t, (22) 

where A = fo(t)\o + /i(i)Ai, and Rq and R\ are defined in the same functional 
form as in Eq. (fP3|) . but in the new basis referred to the crossover point of the total 
non-Hermitian Hamiltonian, Tt T (t). 

The Hamiltonian Ti e f has a complex energy spectrum given by E± = A ± R, 
where 

R = \/g 2 — 2gJ cosa + J 2 , (23) 

and we set J = f (t)R , g — g{t) —i5(t) = fi(t)Ri. From (f2"3"| we obtain the energy 
gap as 

| AS | = 2\ v / g 2 - 2gJcosa + J 2 - 5 2 - 2iS(g - J cos a) \ . (24) 

In Eqs. (|2"3"|) , and below we do not indicate in some expressions the explicit 
dependences on t. 

As can be seen, | AE\ vanishes at the exceptional point, defined in the parameter 
space (S,g) by g 2 + 5 2 = J 2 , g — J cosa. This yields cosa = \Jl — (S/J) 2 . From 
here it follows that the complex energy does not become degenerate during the 
evolution of the system, if S(t) > J(t) for any < t < r. Otherwise, an exceptional 
point appears for some time t < r. 

From Eq. (f2~4")) . taking into account that for many-qubits system sin a ~ 2~™/ 2 
and assuming that S(t)/J(t) 2~™/ 2 for any < t < r , we obtain 

\AE\ > \AE\ min «2min( A /(5- J) 2 + 5 2 ), (25) 

where n is the number of qubits. As can be seen, for any moment of time, < t < r, 
the minimum energy gap |A_B| m i n w 2 min(^/ (g — J) 2 + S 2 . The complex energy 
gap is controlled by the parameter 5(t), and, thus, the non-Hermitian adiabatic 
optimization does not suffer from the typical exponentially small energy gap of the 
Hermitian adiabatic optimization. 
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A rough estimate of the time required for the non-Hermitian quantum optimiza- 
tion can be obtained applying the criterion of Eq. l[2"l"j) to the effective Hamiltonian 
Hef- The estimation of the matrix element \(ip e (t)\H r (t)\ip g (t))\ yields 

i / r , s.„~, / \ i i / \ \ i max I Jo — qj\ sin a . _ 

HMK Ve(* Wr(*) Vff *) » |A F | " ^ 

From here, with help of Eq. (12ip and taking into account that for many qubits 
system sin a w 2~™/ 2 , we obtain 




V. 

T ' 1 



0.2 0.4 0.6 0.8 1 

S 

Fig. 1. Non-Hermitian adiabatic quantum optimization (<5o ^ 0): \AE\/J* as function of the 
dimensionless time s = t/r. The energy gap is rcnormalized and remains non-exponential for any 
time < s < 1 (solid line, So = 0.25J,; dash dotted line, <5o = 0.5J*; dotted line, <5o = J*). 
Hermitian adiabatic quantum optimization (dashed line, <5o = 0): The energy gap between the 
ground state and the first excited state becomes exponentially small at the moment of time 
s = 1/2. 

As an illustrative example, we consider the non-Hermitian adiabatic quan- 
tum optimization algorithm realized by the following linear interpolation: g(s) = 
g(s) — iS(s) = (J* — i#o)(l — s ) and J(s) — J*s, where s — t/r denotes dimen- 
sionless time. Substituting g(s) and J(s) into Eq. (|25p . we find that the energy 
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gap, \AE\, is bounded from below by |A-E| m ; n = 2J r S / y/S^+AJ 2 . At the criti- 
cal point the complex energy has non vanishing gap controlled by the parameter 
So- (See Fig. [TJ) Applying the criterion of validity of the adiabatic approxima- 
tion (|2ip . we find that the evolution time r must satisfy the following condition: 
T > t = 2~"/ 2 J*(5§ + Jl) 1 / 2 /\AE\% lin . From here, in the limit S < J*, we obtain 
t»2-"/ 2 J2/^. 

This inequality means that the time of passing the energy gap point becomes 
non-exponentially small when the number of qubits, n, is big enough. There are two 
main reasons for this: (i) due to the non-Hermitian part of the total Hamiltonian, 
the energy gap is renormalized and becomes non-exponential, and (ii) at the same 
time, the matrix element in (|26p of transition between the ground state and the 
first excited state still remains exponentially small. We should also mention that 
the non-Hermitian part of the Hamiltonian leads to the final life-time of qubits. So, 
the following inequality must be satisfied: 2~™/ 2 J 2 /<5o <C r <C l/<S gu fcit, where S qu bit 
is the charcteristic width of the intrinsic energy level (s) of the individual qubit. We 
should note that the relation between S qu ut and <5o will require a concretization of 
the quantum computer architecture in (|16jl . 



4. Concluding remarks 

In conclusion, we have demonstrated that our adiabatic quantum optimization al- 
gorithm, based on the use of non-Hermitian Hamiltonians, can significantly reduce 
the time needed for optimization of complex combinatorial problems. The main 
ideas of our approach are the following. One uses as an individual qubit a two level 
system which has simultaneously the discrete (intrinsic) states and a continuous 
part of the spectrum. So, the Feshbach projection on the intrinsic states results in 
the finite widths, S qu bit, for of the intrinsic states, which characterize the strength 
of the interaction of a qubit with a continuum. This interaction should not be too 
strong, as it characterizes the life-time of a qubit, ~ ^-/S qu ut- On the other hand, 
the parameter, So, which characterizes the non-Hermitian part of the total Hamil- 
tonian is responsible for the renormalization of the energy gap between the ground 
state and the first excited state. Then, the energy gap becomes non-exponential, 
but the matrix element which is responsible for a transition between the ground 
state and the first excited state still remains exponentially small. All these factors, 
taken together, result in the improving of the adiabatic condition, if the number of 
qubits is big enough. The possible candidates for impleme ntation of this app roach 
could be, for example, the superconducting phase qubits | 20 | 21 | 22 | 23 | 24 | 25 ] w .j ) j c j 1 
are effectively used in recent experiments on quantum computation. In order to 
establish the explicit relation between the parameters 8 qu ut and So one needs to 
build a concrete architecture of the quantum computer model, which in this paper 
is represented in a general form by Eq. (jTSJ) - 
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